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Abstract 

^vq I We analize the renormalization group equations of supersymmetric QCD with A^ = 1 
for the evolution of parton distributions. For this purpose we develope a simple recursive 

%^J . algorithm in x-space to include both regular regions and supersymmetric regions in the 

j~^ I evolution in the step approximation. Supersymmetric distributions are generated within 

CN ■ a radiative model, with vanishing initial conditions for the superpartners. Here we focus 

S^ . on a scenario with broken susy, characterized by a lighter gluino coupled to the standard 

(^ \ evolution and a decoupled scalar quark. Predictions for the all the distributions are 

O ' presented. 
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1 Introduction 



The study of the QCD scahng violations is an important chapter of high energy physics and a 
very important tool for the analysis of future data at the LHC. At large energy and momentum 
transfers, the underlying quark-gluon dynamics is light-cone dominated and controlled by a 
mechanism of collinear radiative emission described by logarithmic corrections to the lowest 
order cross section. In this picture, based on the parton model, initial quark states are assumed 
to be massless and the running of the coupling is linked to the number of flavours nj included 
in the evolution. Crossing intermediate thresholds opens up new channels and new dynamics. 
There is a widely used formalism of perturbative QCD that we think is worth to extend to the 
super symmetric case and which might be useful for experimental searches of supersymmetry at 
the LHC. The picture is particularly appealing if a "supersymmetric content " of the proton 
is found. However, predictions tied to this picture may be used to rule out a possible light 
gluino and/or a light squark. Other applications involve the study of nucleon collisions with 
the atmosphere as in Ultra High Energy Cosmic Rays, where the center of mass energy of the 
primaries can reach several hundreds of TeV's. 

This formalism involves the study of scaling violations induced on the initial state by the 
opening of supersymmetric channels prior to reaching the hard scattering phase. The analysis 
requires the notion of supersymmetric parton distributions, supersymmetric factorization for- 
mulas for supersymmetric DIS and (susy) hadron-hadron collisions on which we elaborate in 
detail. We remark that supersymmetric versions of parton distributions are easy to define, by 
a natural generalization of the usual DIS approach. The phenomenological validity of these 
extensions are clearly linked to the possibility of detecting scaling violations in precision mea- 
surements of final states in p-p collisions, for instance in Drell Yan, or in other processes with 
a distinct final state. With this analysis in perspective we start developing accurate tools that 
we will use in a separate work for a quantification of rates for gluino/squarks effects at the 
LHC. Although there is no evidence of supersymmetry at current energies, or of supersymme- 
try coming from the initial state, quantifying with accuracy these effects we believe, is still an 
interesting task. These studies are also of direct theoretical relevance, since they tell us more 
specifically how to merge the usual parton model dynamics with the new elementary states 
that supersymmetry predicts. 

In this paper we start analizing these issues and make a first step toward quantifying the 
impact of supersymmetry in the initial state. Our analysis, in this paper, is focused on the 
supersymmetric DGLAP evolution, which is more involved compared to ordinary QCD. We 
solve the equations using a recursive algorithm that we formulate and that we have tested 
which allows to perform a direct match between the various regions of the evolution as we 
increase the factorization scale and allow for new (supersymmetric) channels. Part of the 
analysis is of technical nature and several sections deal with the implementation of the method. 
We then study in 2 final sections the implications of the susy evolution and compare it to the 
standard QCD one. Other aspects of the evolution, including some applications for collider 
processes will be analized in a companion paper. 



2 The Method 

Before moving to Supersymmetric QCD (SQCD), we briefly illustrate the method as it applies 
to the case of ordinary QCD. This simpler case will help us establish notations that will be 
later extended to the supersymmetric evolution. We build on previous work of Rossi , which 
has been used before in the previous literature [^, |^, although never fully documented in its 
code implementation. We combine that method, originally suggested for QCD, with recent 
approaches that use an analytic evaluation of the convolution integrals ("weights") [0, and 
generalize it to SQCD. A recently application to QCD in next-to-leading order (NLO) of this 
approach has been implemented by Chuvakin and Smith [0] in their analysis of an evolution 
scheme with a varying number of flavours. To make the discussion self-contained, we have 
elaborated in some detail on the method in an appendix. 

The two-loop running of the coupling constant is defined by 

«(Q^)_2 1 f^ A lnln(QVA^) ^ 1 ,^^ 



27r /3oln(Q7A2)V /3o HQ^A^) 'ln\QyA^] 

where 



Po = yC'c - -^TRHf 



Pi = Y^h-fcGnf-2CFnf, 



where 



N'^ -1 1 

Cg = N, Cp = ^^^, Tn=- (3) 

and N is the number of colours while Uf is the number of fiavours. 
The solution for the running coupling is given by 

ZTT 

with a(Qo) — '^(O)) ^iid Qo denoting the initial scale at which the evolution starts. The 
evolution equations are of the form 



,2 d , „2N «(<5^^ 



^dQ^^'^'''^'^ ^ -^P(-)i^^^iQ))®X^{x,Q'), 



(5) 



with 



Mx, Q') = g.(+)(^, Q') - —q^-'Hx, Q') (6) 



for the non-singlet distributions and 



dr 



(7) 



for the singlet sector. 

We have defined, as usual 

qt^ = qi-q^, qt^=qi + q^, g(+) = s = ^gj+). (8) 

We introduce the evolution variable 

which replaces Q^. The evolution equations are then rewritten in the form 

j^qt\t,x) = (^p(°)(x) + ^/?(_)(x) + ...)®gJ-)(t,x) (10) 

Q'^M^^Q') = (p^'\x) + ^R^+){x)\(^x^{x,Q'), (11) 



In the new variable t, the kernels of the evolution take the form 



i?(±)(x) = Pl^ix) - ^Pl'\x) 

R{x) = P«(x)-Ap(o)(^). (13) 



(-) 



Equations (0) and ( [Tl| ) are solved independently for the variables g) and Xi respectively. 
Finally, the solution q^~^^ of eq. (|12D (or the singlet equation) is substitued into Xi iii order to 
obtain g| . 

The equations can be written down in terms of two singlet evolution operators E±{t, x) and 
initial conditions q±{x,t = 0) = q±{x) as 

^E± = P±®P±, (14) 

whose solutions are given by 

qt\t,x) = E(^_)0qt^ 

Xiit,x) = E(^+)(g)Xiix). (15) 

The singlet evolution for the matrix operator E{x,t) 



EGF{x,t) EGG{x,t) 



(16) 



^^^ = P^E{x,t) (17) 

is solved similarly as 



S')^-''-'«(':;;i;' 



The unpolarized leading order kernels are expanded in a as 

P^,ix,as) = (^) P,f(x) + {^yPS\x) + - (19) 

Their LO expressions are given by 

for the non-singlet sector, and by 



^qq l-^J ~ ^qq,NS 

PJ^^g^{x)=2TRnf(x^ + {l-x 
l + fl-a;)2 
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P^q\-)=C 

P(;)(x) = 2iV, ((^^ + ^ - 2 + x(l - x)) + 1^(1 - x) (21) 

in the singlet sector. Some simple identities for the "plus' (+) distributions and their definition 
are given in the appendix. 

In the actual numerical solution of the equation, one would like to have at hand a recursion 
relation which can be implemented in a computer program dynamically at run-time. In order 
to develope recursion relations we proceed as follow. We first observe that as far as we are not 
resumming double logarithms in Q and x, the solution eq. (18) can be expanded as a series of 
convolution products 



q{x, t) = qo{x) + tP ® qo{x) + —P ® P ® qo{x) + ... (22) 



and transformed into a recursion relation for some coefficients Cn{x) 



Co{x) 

Cn+1 



6{x-l) 
P ® Cn{x) 



equivalent to the expansion for the evolution operator 



E{x,t) = Y.CN{x) 



n=0 



t" 



n\ 



(23) 



(24) 



Written as a recursion relation, the equation is as easy to implement as a calculation in moment 
space, but with no need to perform any moment inversion. The possibility of using this expan- 
sion as an alternative way to evolve parton distributions fast and with accuracy especially for 
larger sets of coupled equations, such as for supersymmetric theories, is the subject of the follow- 
ing sections. Specifically, Rossi's ansatz P], originally introduced in the context of the photon 
structure funtion, differs from (^) only by an over all normalization Cn = (~1)^ (2//3o) 



\N 



A^l 



:^) 
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n=0 



n\ 






(25) 



Connecting this expansion to other expansions is also pretty straightforward. A very elegant 
basis in which to expand is the Laguerre basis, introduced by Furmanski and Petronzio 0]. A 
relation between Rossi's basis and the FP basis can also be derived but the result is not 
particularly illuminating. 

Inserting the expansion ( ^Sf ) into the evolution equation one gets some recursion relations 
for the functions An+i{x) in terms of the An{x). These are obtained by comparing left hand 
side and right hand side of the evolution equations after equating the logarithmic powers with 
a running strong coupling constant a((5^)- A pretty detailed study in the polarized case can 
be found in ref. [^j. For computational purposes, the recursion relations for the evaluation of 
the An{x)^s can be implemented in various ways. The one that we have implemented involves 
the direct solution of the recursion relations as illustrated below. In leading order in a(Q^) in 
QCD we get for the unpolarized kernels 



A 



v,± 

n+l 
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int the non-singlet case and 
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ln(l 



X] 



Po Jx y 



X y y Po 

dyyAl+{y)-xAl+{x) 



y-x 



1 dy ] x {x-y) 

X y \ y y 



+ iM^(y) 



(27) 



^^+i(^) = ^C, I ^ ^"^:, "" At{y) 



/3o Jx y z 



(3o Jx y y-x (3q 

- ^r^|i-2+.(i-.)U„^(y)-A„«(x) 

Po Jx y iz J 

in the singlet. Extensions to the NLO case of this approach are pretty straightforward, and the 
procedure will be more clear once we will discuss its implementation to solve the Supersym- 
metric DGLAP (SDGLAP) equations. 

3 Supersymmetric scaling violations 

In A^ = 1 QCD gluons have partners called gluinos (here denoted by A) and left- and right- 
handed quarks have complex scalar partners (squarks) which we denote as qi and qn with 
q = qi + (JR (for left-handed and right-handed squarks respectively). 

The interaction between the elementary fields are described by the SU{3) color gauge in- 
variant and supersymmetric lagrangean 



+q,I Pq, + D^qnD^f + D^qLOf^f + igV2 (^X-nqlT^qu + X'^LqlnT'^qm - h. c.) 



'29^[(lLiT^(fLi-qRiT''q\i,j +mass terms, 



(28) 



where a runs over the adjoint of the color group and i denotes the number of flavours over 
which we sum. 

Past studies of these effects lead to the conclusion that the information that supersymmetric 
initial states carry along the evolution can be easily absorbed into scaling violations coming from 
ordinary QCD. These studies require the knowledge of the matrix of the anomalous dimensions 
(for all the moments), or of the corresponding Altarelli-Parisi (DGLAP) kernels. We recall that 
the leading and next-to-leading anomalous dimensions are known p[ since long ago, at least in 
the case of partial supersymmetry breaking, in a scenario characterized by a light gluino and 
a decoupled scalar quark (we present in an appendix the leading order form of these kernels 
for completeness). In this work we consider the possibility of a radiatively generated gluino 
distributions, in analogy with the case of standard QCD for the gluons. In fact, in QCD one 
can use a simple model of the distributions at low Q {Q = Qo) and run the DGLAP equations 
up in energy in order to generate distributions of gluons at any higher scale Qj. Of course, it 
is well known that in QCD such initial evolution scale is model dependent. Various amplitudes 
at higher energy emerge, depending on the underlying assumptions on the form of their initial 
shapes. 



This approach has been widely used in the hterature and can be a way to connect low 
energy quark models -which have no gluons but a phenomenological confining potential - to 
true QCD scattering amplitudes. One can assume a zero distribution of glue at the lowest scale, 
or a model dependent non zero one, and then use the evolution equations to dress the matrix 
elements describing the distributions by logarithmic corrections. Once the final scale -here 
identified as the factorization scale of the process - is reached, the distributions are convoluted 
with the usual -on shell- parton cross sections to generate the full hadronic cross section. 

In hadronic collisions, Qf is usually a fraction of the center-of-mass energy, or a fraction 
of the large pt of the final state jets. It is not uniquely identified and a different choice of Qf 
underscores a drastic sensitivity of the expansion on this scale. 

Alternative choices for density of the constituents at low Q are also possible, but, ultimately, 
whatever the choice for the parton structure, it has to match the scattering data available from 
DIS and colliders experiments. In the case of exact SQCD, it is natural to parametrize the 
parton distributions, now with a susy content in the form 

q^Lix,Q^) + l,Lix,Q')), (29) 



where the new elementary constituents (squarks and gluons) can be vanishing at low Q and be 
radiatively generated by the evolution, similarly to the gluon case (in standard QCD) contem- 
plated above. In the case of a broken susy, the form of F2 does not change from the standard 
QCD form. 

A gluino or a squark parton distribution (a more detailed analysis will be presented else- 
where) is the exact correspondent of a gluon or a quark parton distribution, namely a light- 
cone dominated diagonal (non local) correlation function in spacetime, Fourier transformed to 
(Bjorken) x-space. 

Similarly to the QCD case, in the case of exact A^ = 1 supersymmetry we define singlet and 
non-singlet distributions 



qv{x,Q^) = ^(gi(a;,g2) -gi(x,Q^ 
1=1 

qv{x,Q^) = ^(gi(a:,(5^) -^i(a;,g^ 
1=1 

q+{x,Q^) = Y.[q,{x,Q^) + q,{x,Q^ 
1=1 

q+{x,Q') = Y.i'i^i^^Q') + Ux,Q')). (30) 



1=1 



The evolution equations can be separated in two non-singlet sectors and a singlet one. The 
non-singlet are 






Q^-J7^QV{X,Q^) = -^^{Pqq®qV+Pqq®qv) 






{Pqq ®qV + Pqq O Qv) , 



(31) 



and the singlet, which mix qy and qy with the gluons and the gluinos are 



Q' 



dQ^ 



X{x,Q') 
q^{x,Q') 
q+{x,Q') 



Pgg 


Pgx 


Pcq 


Paq] 




r G{x,Q^] 


PXG 


Pxx 


Pm 


Px^ 


® 


\{x,Q^) 


PqG 


P,x 


Pqq 


Pqs 


q+{x,Q' 


PsG 


PsX 


Pqq 


Pqq . 




[t{x,Q' 



(32) 



There are simple ways to calculate the kernel of the SDGLAP evolution by a simple extension 
of the usual methods. The changes are primarily due to color factors. There are also some 
basic supersymmetric relations which have to be satisfied that will be analized below. They 
are generally broken in the case of decoupling. We recall that the supersymmetric version of 
the (3 function is given at two-loop level by 



/5f 



i (IICa - 2nf - 2nx) 

- (SACl - lOCAUf - IOCaUx - eCpUF - GCxnx 



(33) 



with Cx = Ca = Nc for the case of Majorana gluinos and the ordinary running of the coupling 
is replaced by its supersymmetric running 



c^HQl 



2n /?ifln(gVA2 



1 - 



/5f lnln(gVA2 
'P! In(QVA^) 



+ 0{ 



ln'(QVA2) 



(34) 



The kernels are modified both in their coupling {a -^ a'^) and in their internal structure 
(Casimirs, color factors, etc.) when moving from the QCD case to the SQCD case. In our 
conventions an index "S" stands for a supersymmetric component (regular, i.e. non supersym- 
metric, kernels do not carry such an index), but we will omit it when obvious. 



4 Evolution and Matching 

Susy is necessarily broken in the real world and therefore, the way the breaking occurs dictates 
both the mass spectrum and guides the structure of the QCD scaling violations as well. There 
are several parameters that appear in the evolution, mx and rriq, the masses of the gluino and 
squarks, hence a complete analysis includes various scenarios, on which we briefly elaborate. 
In a realistic scenario with a broken susy, the squarks have much larger mass compared to the 
quarks and the gluinos have a Dirac or Majorana mass mx (or a combination them). In our 
case, in order to establish the evolution scales which are of phenomenological interest, it is 
convenient to assume 1) that the scalar quark decouples from the evolution and 2) that the 
2-gluino production threshold ■m2x = 2mx is the intermediate scale, separating in Q'^ the regular 
QCD region from the supersymmetric one. In this scenario some of the splitting functions - 
specifically Pqx and Pxq - are zero, since no collinear emission is associated with massive partons, 
unless the symmetry is effectively restored. 
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We impose the separation condition 

Qo < m2x < Q 
and only one non singlet evolution equation is considered 



(35) 



r>2 d , ^2^ a(Q^) 



Pqq®qV- 



(36) 



The simplified singlet equations, which mix qv and qv with the gluons and the gluinos 
become 



Q' 



d 



dQ' 



A(x,g2) 
q{x,Q^) 



'gg Pgx 


Pcq] 




\G{x,Q' 


AG -PaA 


Px, 


® 


A(x,g2 


'<jG PqX 


p 




[ q{x,Q': 



Formally, we can write down the expression of the complete solution in the form 



qNs{x,Q}) = q{x,Ql)+ dlnQ^ PNs{x,a{Q'^)^qNs{x,Q^ 



+ 



'^Ix 



dlnQ^ PNs{x,a{Q^))(g)qNs{x,Q^) 



for the non-singlet equation and 



Gix,Q}) 
\ix,Q}) 
q{x, Q}) 



Gix,Ql] 


qix,Ql) 



+ 



d\nQ''P{x,a{Q))0 



Gix,Q^) 



g(x,g2) 



+ 



1^' d\nQ^P^{x,a^{Q^))® 



G(x,g2) 



g(x,g2) 



(37) 



(3J 



(39) 



for the singlet one. In analogy with eq. (|25|), now we introduce supersymmetric coefficients 
A"'{x) beside the usual non supersymmetric An{x), and impose recursion relations on the initial 
ansatz q{x, Ql) of the form 



A^ix) 

aNS 



A^ix) 



S{l-x)0qix,Ql) = q{x,Ql), 

~^Pqq yS) An, 

Po 



n=0 



n! 



a(go) 



A, 



NS 
n+1 



— P^ 6d A^^ 



(40) 
(41) 

(42) 



(43) 
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n = 1, 2, ..., which solve the equation in the non-singlet sector. 
In the singlet sector we get the solution 



G{x,Ql) 



G{x,Ql 


q{x,Ql\ 



OO -| 

+ E 



n=l 



n\ 



+ 



OO -| 



n=l 



n\ 






log" 



a{Q) 



a[m2x 



constructed recursively through the relations 



A'oix) 
A',ix) 
Alix) 

An+l{x, Q ) 



Gix,Ql) 


qix,Ql) 



A" 



n+l 



X] 




P® 



A'oi 



x] 



A',ix) 
Alix) 

^n+l\X) 
An+l{x) 



00 -| 
n=0 "■• 






X 



XK' 



a4„(x) 
log" 



A^^ix) 



. a(Qo) . 



log" 



' a(m2A) ' 
. a(Qo) . 



(44) 



(45) 
(46) 

(47) 



(4J 



77. = 1, 2, ..., having used a vanishing gluon density at the starting point of the supersymmetric 
evolution. 



5 Supersymmetric Relations 

An exact supersymmetric scenario is probably interesting only for analizing theoretical issues 
concerning the evolution, or to study the shape of the distributions at extremely high energies, 
when, again, all the supersymmetric partners effectively become mass degenerate. This scenario 
can take place at few TeV's or at several TeV's, depending upon the assumptions underlying 
the way supersymmetry is restored. Here we simply focus our analysis on a light, up to an 
intermediate-mass Majorana gluino. Other aspects of the evolution, such as the interplay 
of Dirac and Majorana gluinos and the impact of various patterns of susy breaking will be 
considered elsewhere. 

As we have already mentioned, the study of the anomalous dimensions - in leading order - 
for A^ = 1 QCD can be found in |l[] and now we are going to elaborate on that. 
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We recall that in a regularization scheme which is manifestly supersymmetric there are some 
supersymmetric relations which are satisfied by the anomalous dimensions. This is true, for 
instance, in the Dimensional Reduction DR scheme. We recall that to leading order the MS 
scheme and the supersymmetric DR scheme give coincident results. We recall that in the DR 
scheme the traces are kept in 4 dimensions, and the chiral projectors are treated as usual, with 
a completely anticommuting 75. The loop momenta are evaluated in n-dimensions. In other 
schemes, such as the t'Hooft-Veltman scheme, 75 is instead partially anticommuting. The MS 
scheme for chiral states is usually based on this second definition. The relation between the 
two schemes is to leading order 

MS DR ^^> 

valid both for the non singlet and the singlet anomalous dimensions. 

The same is true also for the factorization scheme dependence of the coefficient functions. 
As for the coupling constant, the definition of a in the two schemes, the expression of ajjg can 
be used also in the DR scheme as far as the Aqcd scales of the two schemes are related by 
A;^^ = Ajjgexp (Ca/QPo)- This last change is tiny and will be neglected. 

Part of the supersymmetric kernels can be obtained at this order by a simple change of 
group factors, from the standard QCD kernels. For instance, the substitutions Cp -^ C\ 
{Cp = (A^c ~ l)/(2A''c)) and n/ -^ n\, allow us to obtain the kernels {\\, Xg,g\) from the 
ordinary kernels {qq,qg,gq)- We use C\ = Ca = Nc, the number of colors. As for the choice 
of the type of representation for the gluino ( Majorana or Dirac), we recall that in the Dirac 
case we set nx = 2Ca while in the Majorana case nx = Ca- The organization of the terms in 
the splitting functions may differ from reference to reference, due to the various manipulations 
one can perform on the plus distributions. 

The supersymmetry relations are given by 

(50) 
and by 



P99 + ^Ag ■ 


— Pg\ + P\X 




= Pgs + Pxs 


gg ' sg 


= PqX + PsX 


^qq ' ^sq 


-* <js ~r J ss 



(51) 

In general, for a decoupled scalar quark, the only symmetry one would expect is eq. (^Up. In 
the case of Majorana gluinos, for a scenario with a decoupled squark, it is interesting to observe 
that this relation remains valid for x < 1 as well. It can be extrapolated to include the x = 1 
point in the case of zero number of fiavours (supersymmetric gluondynamics) . 

The evolution has to respect - both in the case of exact susy and of susy breaking - 1) baryon 
number conservation and 2) momentum conservation. There are two sum rules associated 
with these conserved quantities, which are generally used to constrain the phenomenological 
parametrizations of the distributions in direct applications. Below the m2x threshold we need 
to satisfy the usual QCD relations 

/ dx (xG{x, Q2) + xq'^+) (x, Q2)^ _ I (52) 
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for momentum conservation and 

/ dxq^'\x,Q'^) =3 (53) 

Jo 

for baryon number conservation. These require that the anomalous dimensions satisfy the 
relations 



/ XdX {Pgi{x) + Pgi{x)) = 

Jo 

/ dxPNs{x) = 0, 

JO 



(54) 
respectively, where i = q,g. In leading order the second equation is simply 

/ dxP^f{x) = 0. (55) 

Jo 

Moving above the 2-gluino threshold the momentum sum rule becomes 

dx {xG{x) + x\{x) + xq'^'{x)\ = 1 (56) 



for momentum conservation, while eq. (|53[) remains unaltered. We get the new momentum sum 
rule (or second moment sum rule) 

;.<;.(P.M + n.M + a.) = 0, (57) 

with i = q,g, X. The supersymmetric version of eq. (|55D is simply obtained by replacing Pjg^(x) 
by it supersymmetric counterpart P^J-^^x). It can be checked easily that the kernels in the 
appendix satisy these relations. 

In the case of exact supersymmetry we need to keep into account the scalar quark contri- 
bution in both equations. In particular, the equation for the second moment becomes 



"1 
/o 
which implies that 



xdx (xG{x) + xX{x) + xq^~^'{x) + xq^~^'{x)) = 1, (5J 



/ xdx {Pgi{x) + Pqi{x) + Pxi + Pcii) = 0, (59) 

Jo 

with i = q,g, A, q. As for baryon number conservation, equation (|53|) gets modified into 

(60) 





/' 

JO 


dx (q^ 


\x)+q^- 


-\x)) - 


= 3 


and 


using (|31) one gets two relations 




JO 


d^ K 


+ Psl) 


= 






JO 


dx (Pf, 


+P.I) 


= 





(61) 
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(below we will drop the susy index S in front of the kernels when obvious). The first relation 
in the equation above, for instance, clearly implies that the end point contributions in the 
ordinary Pgg kernel are to be modified in order to insure conservation of baryon number. In 
the case of a susy breaking scenario such additional end-point contribution is trivially absent 
and the form of the Pgg kernel remains the same, except for the replacement of a with a^, the 
super symmetric coupling and Pq with /3^. 

6 Applications 

In order to apply the formalism to SQCD we proceed by discussing the method for a scenario 
with a broken susy (decoupled squark). We proceed from the non singlet case, analizing in 
details the intermediate steps of the algorithm. 

We start with a solution of the standard non-singlet DGLAP equation assuming the bound- 
ary condition Aq{x) = q^{x, Ql) and run the equation up to the gluino threshold and construct 
the coefficients recursively. We arrest the coefficient up to a desired order (n), which can be as 
large as 30. In general, the rate of convergence of the asymptotic expansion changes with the 
value of the momentum Q. 

The solution is then constructed in the region Qo < Q < m2\ as 

q x, m2x) = 2^ j— In 62 

^0 ^' V «Wo) / 

and used as initial condition for the next stage of the evolution, which involves the region 
'^2A < Q < Qf, with Qf being the final evolution scale. 

At the next stage, we set Aq{x) = S{1 — x) ^ q^{x,m2x) and solve recursively using the 
super symmetric version of the kernels. The strong coupling constant Q;(<5^) and its running 
are replaced by their supersymmetric version a^lQ"^), and so are the coefficients of the beta 
function (/3j -^ f3'-). Finally the solution is written down in the form 

,.(,,Q^) = ±^io,"(^) + i^iog'(^) (63) 

^0 n! V "(Qo) J ^i n\ \ a{Q) J 

where n' is the index at which we arrest the supersymmetric recursion. In our implementa- 
tion we have kept the values of n and n' very close. 

In the singlet case the procedure is not much different. We evolve according to the standard 
DGLAP equation in the region below the 2-gluino threshold, having set to zero the contribution 
from the gluino at the beginning (up to the supersymmetric threshold). This is equivalent to 
having set A^[x) = for any n, which means that in the region below m2\ there is no radiative 
production of gluinos in the initial stage. We iterate the recursion relations up to a given value 
n of the index n with the standard DGLAP. The initial conditions for the next stage of the 
evolution are then fixed by the relations 

Al^{x) = g+(x,m2A), 

= i:^i°g"(^). (M) 

^0 "-' V "Wo) / 

Af^{x) = G{x,m2x) 
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Ki.x) J „ ( a{m2x) \ 



A',{x) = X{x,Ql) = 0. 

(66) 

After this, we determine recursively the coefficients An of the supersymmetric expansion 
,,+ _ ^Cf /I dy yAl+iy) - xA^+ix) 

^n+l — 



j3§ Jx y y-x 

'^-. log(l - x)AlHx) + ^ /' ^(1 + z)Al-iy) - ^ir(x) 



Po Po "'^ y Po 



riv(i-2^ + 2-)<fe). 



/^o -^^ y 



^x . ^ _ .C'a f^dy yA^iy)-xA^^{x) 



2 ^ f^dy [2 ^ \ 7„, ^ 2 ^ f^dy [2 

' Po Jx y y-x Pi 

-4^ /' ? (7 - 2 + ^(1 - ^)) i^„(l/) - i^Jx). 
Pq Jx y \z / 



Finally, we construct the solution in the form 

;(,. Q^) ^ I: ^ log" f ^) + i ^ log" f ^) 



(67) 



Po "'^ 1/ Po 

-|n. /' ^(1 - 2z + 2z'))AUy) - 4§i^(x) log(l - x), (6J 

Po "'^ 1/ Po 



(69) 



(70) 



(71) 



where we have arrested the supersymmetric recursion up to the index n'. Here f{x,Q ) indi- 
cates a singlet quark, a gluon, or a gluino distribution, with A(^{x) and A{(x) denoting their 
corresponding coefficients in the expansion. Singularities emerging from the lower integration 
point {y = x) are the tricky part of the game, as expected, but can be handled with various 
techniques. They will be discussed briefly in the last section. 
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7 NLO Extensions 

Let's now move to a next-to-leading order (NLO) analysis of the evolution. Our discussion here 
is partial and does not include contributions due to the emergence of new anomalous dimensions 
as we move across the supersymmetric threshold. In recent work it has been shown that in 
the ordinary distributions of quarks and gluons of QCD these effects are important, especially 
at small-x. These changes are expected to produce only a slight modification of the algorithm 
presented below, and simply amount to a modification of the boundary condition as we move 
across the m2x point. They will not be analized any further in this work and will be assumed 
to be negligible. As a second point, we remark that the extension of the procedure outlined 
below is easy to generahze to the more general case of exact susy. 
To NLO the ansatz becomes 

and inserting the usual running of the coupling 

we get the recursion relations 

2 
An+i = -—An{x) 



A 







5„+i(x) = -Bn{x) - ( ^A„+i(x) 1 - j^P^'^ ® A„(x) - ^P(°) ® 5„(x) 

V^Po / 47rPo Po 

(74) 

which are solved with the initial condition Bq{x) = 0. The initial condition for the An{x) 
coefficients (i.e. Aq{x)), is specified as in the previous section, with q{x,Ql) identified as the 
leading order ansatz for the initial distribution, i.e. 

Ao{x) = 6{1 - x) ® q{x, Ql) = g(^°)(a;, Ql) (75) 

Running the RGE's below the gluino threshold region (up to Q = 'm2x) and solving the 
recursion relations (|74D, we get the solution (arrested at a recursive index n) 

^0 ^- \a{Qo)J ^0 n\ \a{Qo)J 

which is used to fix the intial condition for the second stage of the evolution, the supersymmetric 



one 



q{x,m2\) = q {x,m2\) + a{m2\) q {x,m2x), (77) 
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with 



LO( ^ ^AJ^ (a{m2x) 

q {x,m2x) = 2^ j— log 



n\ "° \ «((5o^ ' ' 



„=o '*• \ "WOJ 

n o 
^ -Dri 

n=0 '^^ '"'' V "('^o) J ' 

(78) 



g^^^(x,m,,) = J2^iog-(^^^] 



The supersymmetric recursion relations are then given by 



io(x) = 6{l-x)^q^^{x,m2x), 

Bo{x) = g^^^(x,m2A), 
2 



i„+i(x) = --^P(°)^®i„(x) 



^n+l(x) = -5„(X) - (^^+l(^) j - ^^^'^" ® ^(^) - ^P^'^' ® ^n(x) 



-^n(x) + ( :^P^'^' ® in(a:)) - ^P(^)" ® in(x) - ^P^'^' ® 5.(^) 



A' 

.2/55 



(79) 



We finally construct the general solution in the form 



,(x. Q^) = ,(x. ,„«) + t M£) ,„g.. (^) + a(Q) E ^ log" (^) . (80) 

^0 n! V«("^2a)/ ^0 ^' \a{m2x)J 

Eq. (^) can be expanded in components, since it is valid in matrix form 

2Po 



2 



4/3o^ 



and similarly for the evolution of the gluon density 

^^+i(^) = -mx) + ^{p!^'J'®Al^x) + P!^'J'®Ai{x) + P^'^'^A'^{x: 

2Po 

Po 



f82) 
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while the gluino density is obtained using the recursion relations 

?5 






2 



^ wr ® ^r(x) + pff ® s^(x) + Pi?" ® p„\x)) . 



(83) 



Notice that as initial condition for the gluino distributions we take an identically vanishing 
function at the scale Q = m2\ both in leading and in next-to-leading order. 



Ao{x) = 
Bo{x) = 0. 

(84) 

8 Numerical Results 

It is expected that a large gluino mass, for a fixed factorization scale Qf in the evolution, 
lowers the size of the scaling violations and their impact on the supersymmetric cross section. 
On the other side, scaling violations induced by the susy evolution should grow as we raise 
the final evolution scale. Therefore it seems natural to study the effect of the susy evolution 
in two different setups 1) for fixed mx and a varying Qf or 2) for a varying mx at a given 
factorization scale Qf. We have performed both studies and the results are shown in figs 1-14. 
The implementation of the unpolarized first stage (QCD) evolution is performed in the MS 
scheme, which is by now standard in most of the high energy physics applications. We introduce 
valence quark distributions qv{x,Qo) and gluon distributions G{x,Ql) at the input scale Qq, 
taken from the CTEQ3M parametrization |10 



q{x) = AoX^^(l-x)^2(l + /l3X^^). (85) 

Specifically 

xuvix) = 1.37a;°-^9^(l-a;)=^-^^[l + 6.25x°-H 
xdvix) = 0.801x°-^^^(l-x)^-i9[l + 1.69x°-37^] 
xG{x) = 0.738x-°-2^^(l-a;)^-3^[l + 7.30x] 

^u{x)^^d{x)_ ^ o.547a;-°-286(i_a;)«-34[i + i7.5x] 

u{x) + d{x) 
xs{x) = 0.5 a; 

x{d-u) = 0.75 2;^-^^(l-2;)^-3^(l + 30.0 x) (86) 

and a vanishing anti-strange contribution at the input. In figs. 1-13, where we have studied the 
valence quarks, the gluon and the gluino distribution for a varying gluino mass ranging from a 
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light to an intermediate gluino (10-40 GeV) up to a value of 250 GeV. These results generally 
point toward small scaling violations, which become more appreciable as we move closer the 
smaller-x region (in particular for gluons and gluinos). We have chosen as initial evolving scale 
Qi = 2 GeV in all the runs. 

Fig. 1 shows the shape of the distributions close to their initial evolution scale, with Qi = 2 
GeV and Qf = 5 GeV. In fig. 2 we show the scaling violations induced by the susy evolution 
on the u and d quark distributions, from the initial scale up to the final scale of 100 GeV for a 
gluino mass of 10 GeV. In this figure we show the regular versus the supersymmetric evolution 
for these two distributions. The modifications appear to be quite sizeable. The presence of 
supersymmetry in the evolution shows up as a lowering of the maxima of the distributions with 
a shift toward the small-x region. 

A better distinction between the non susy from the susy result is illustrated in fig. 3, 
which shows a comparison between regular and susy evolution for gluons. As expected, these 
differences get more pronounced moving toward the region of smaller x, due to the rise of the 
gluon distribution and to the small-x structure of the kernels. In fig. 4 we show the singlet 
q^^^ distribution. The difference between the regular and the supersymmetric evolution is slim, 
given the low factorization scale (100 GeV) chosen in this run. 

Gluino and gluon distributions are shown in fig. 5. The gluino distribution is approximately 
2 orders of magnitude smaller compared to the the supersymmetric gluon distribution and 
grows fast at small-x. Fig. 6 shows a plot of the gluon distribution for a large final factorization 
scale Qf = 1 TeV and a varying gluino mass (100,150 and 200 GeV respectively). Shown is 
also the regular evolution of the gluon density, which is lower than the susy ones at larger x 
vales and faster growing at smaller x values. Fig. 7 shows the factorization scale dependence 
of the gluino distribution for a sizeable gluino mass (100 GeV) and large factorization scales 
Qf = 1,2 and 5 TeV. We plot in fig. 8 the gluon and the gluino distributions for very large 
evolution scales for a realistic gluino mass of 250 GeV and varying factorization scales. We 
have chosen Qf = 5 and 10 TeV respectively. The dependence on the final scale appears to 
be quite mild. Fig. 9 shows the variation of the gluino distribution for different gluino masses 
(40,100,200 and 250 GeV) and a fixed final evolution scale Qf = 1 TeV. There is a reduction of 
the small-x growth of this distribution at smaller-x values as the mass of the gluino is raised. In 
fig. 10 we show the dependence of all the quark distributions on the factorization scale Qf = 5 
and 10 TeV in the case of a supersymmetric scale m2A of 250 GeV. 

We study the impact of these corrections on future collider experiments by showing results 
for the 2-gluino production in a p-p collision. As can be seen from fig. 11, the production cross 
section is quite small -compared to standard QCD rates-, but gets enhanced by the inclusion 
of susy scaling violations especially at larger energies. We have set the gluino mass at 250 
GeV. In the same figure we show the dependence of the cross section on the factorization scale 
(510,600 and 700 GeV). The dependence is sizeable, although the total rates remain small 
compared to the QCD background at these energies ( 650 GeV - 2 TeV). We have shown in 
figs. 12 and 13 the partial contributions to fig. 11 coming from the g g ^ XX channel and the 
qq channel. The gluon contribution is the dominant part in both the regular QCD case and in 
the supersymmetric case. 
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9 Conclusions 

We have solved the supersymmetric DGLAP equation in a scenario with a coupled gluino and 
a decoupled squark using an algorithm based on x-space. In particular, we have illustrated the 
evolution of the distributions of quarks and gluons and their supersymmetric versions using a 
radiative model. Although the window on a light gluino is now rapidly closing and the hope 
to detect supersymmetry from scaling violations in the initial state at the LHC with a light 
gluino is slim, the possibility of analizing experimentally the impact of heavier supersymmetric 
particles in the initial state remains, however, an important issue. In a specific example (p-p 
— > AA), we have seen that the total cross section is sensitive to initial state susy scaling violations 
and on the factorization scale chosen. For 2-gluino production, for instance, the rate is small, 
much below the usual QCD background, but gets sizeably enhanced when susy is included. 
Would be interesting to see how much we can rise the gluino mass and still obtain a signal on a 
final state which can not be compensated by the usual re-adjustement of the several parameters 
that describe the ordinary QCD parton distributions, or is comparable to it. In particular, the 
strong factorization scale dependence may be reduced by the inclusion of radiative corrections 
in the initial state. This and other related issues will be analized elsewhere. 
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9.1 Appendix 1: The Kernels 

Various relations among different types of "+" distributions can be derived. The kernels given 
in the literature may differ in their final expressions due to rearrangements of the correspond- 
ing "+" functions. We give here the expression of these kernels and illustrate some of the 
manipulations needed to reorganize them in a standard form. 

In the derivation of the recursion relations we have used the identity 

"""< ' V4,fe)=/'*l^i:M^£ii*i + iog(i-.M„w (87) 



Ix y V(l -x/y)+J Jx y y-x 

Similarly, it is not hard to show the identity of the convolution products 

^'^"'^ -Mx) = /' ^ y" + "; AM - A(x) r d.i^ - A(x) '•' '^ ^" + "' 



xl Jx y{y-x)y Jo 1 - z Jx y y'^{y - x) 



with 



'—-] -l-a; + ^<5(l-x) 
x)^ / 2 



A = [^ dy^ yAnjy) - xA^jx) 



X y y-x 



2 log(l - x)Ar,ix) + ^A{x) - f'^(l + z) A,{y) 
2 Jx y 

(89) 
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After some manipulations one can show that the two expressions given above are equal and 
therefore 



l + x^ 



x). 



1 - X + -6{1 - x) 



(90) 



In leading order, the supersymmetric expression of the standard (qq,qg,gq,gg) QCD kernels 



are obtained by replacing a 



a 



and the dependence on the coefficient of the /? function by 



its supersymmetric version Pq -^ Pq. 

The leading order expression of the kernels in the case of a decoupled scalar quark are given 
below. We have omitted the superscrit "S" for simplicity. We remark that the usual QCD 
kernels, after the embedding in the supersymmetric evolution, do not acquire new terms at 
X = 1, except for Pgg. 
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XI 
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Cx 
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99 



Cx 







Pi^ = -f 

nf{l-x) 
Cf 



1 



(1 — x)+ X 



-2 + x(l -x) 



1 -x + -5{l -x) 



x" + (l 
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,(1-^)- 



l-x+-<5(l 



X 
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C. 



1 +x2 



X 
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10 Appendix 2. Discretizations 

For the calculation of the weights used in the numerical analysis we follow closely ref. [Q 
also implemented in 0]. At each x- value, these authors use an approximation characterized 
by weights which are calculated analytically, together with an interpolation formula for the 
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integration function. The method allows to monitor the singularity appearing in the small-x 
region and to achieve a very good numerical accuracy. The method speeds up in time the 
calculation by a large factor, but becomes tedious when moving to a higher order, since all the 
integrals have to be exactly discretized and the logarithms extracted in each sub-interval. Here 
we briefly illustrate the method as it applies to our case. 

We briefly recall the numerical strategy employed in this analysis. We define P{x) = xP{x) 
and A{x) = xA{x). We also define the convolution product 

The integration interval in y at any fixed x- value is partitioned in an array of increasing points 
ordered from left to right {xo,xi,X2, ..-jXnjXn+i) with xq = x and Xn+i = 1 being the upper 
edge of the integration region. One constructs a rescaled array {x,x/xn, ■■■,x/x2,x/xi, 1). We 
define Si = x/xi, and s„+i = x < Sn < Sn-i < •••Si < sq = 1. We get 

,(.).!: r"^(£)p(£)i(,) (93) 

^A, y \yj \yj 
At this point we introduce the linear interpolation 

A(y) = (i- y-^' ) Aix,) + y~^' A(x,+o (94) 

and perform the integration on each subinterval with a change of variable y— > x/y and replace 
the integral J{x) with its discrete approximation Jn{x) to get 



Jn(x) 



A{x,)-^ f '^P{y){y - s,) 
1 - si J Si y 

EAix,)^^ r ^P{y){y - s,^,) 
A—^ Si — SiMi Jsi+^ y 



~i ' ' Si ~ Sj+i Jsi+i y 

- J:A{x,)^^r"^-^P{y){y-s,_,). 
"—{ Si_i - Si Js, y 

(95) 
Introducing the coefficients W{x,x) and W{xi,x), the integral is cast in the form 

n 

Jj^(x) = W{x, x)A{x) + J2 ^i^u x)A{xi) (96) 

1=1 

where 

W{x,x) = -^(''^{y-s{)P{y), 
1 - Si J SI y 

W{x,,x) = -^^r ^ (y - s,^,) P{y) 

Si ~ Sj+1 JSi+i y 

''i-i dv 

\y-s,_,)P{y). 



Si-i — Si J Si y 

(97) 



We recall that 



',.. ^w 



/o (1 - x) 

and that 
1 



1-y 



'l-x) 



Jx y y-x 



as can can be shown quite straightforwardly. 
We also introduce the expressions 
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(98) 



(99) 



Ino{x) 

Jni{x) 
Jnti (x) 
Jnti(x) 



Sl 



1 



■Si — Sj+i 

1 

"Sj — 1 Si 

1 



log(si) + log(l -Si^ 



1 - Si+i \ I- Si S^+l 

log —, + Si+1 log 



1 — s,- 



1 ~ -Sj+l Si 



log I I + Si_i log I -^- I + S,_i 



1-s 



j-1. 



Si~l. 



1 - S,; 



1-S 



logsi. 



Using the linear interpolation formula (Q) we get the relation 



i = 2,3,..A^ 



(100) 



dyyAn{y) -xAn{x) 



X y 



y-x 



- log(l - x)An{x) + A„(x)Jno(a;) 



N 



+ ^An{xi) {Jni{x) - Jnti{x)) 
1=1 



(101) 



which has been used for a fast and accurate numerical implementation of the recursion 
relations. 



References 

[1] C. Kounnas and D.A. Ross, Nucl. Phys. B214:317, 1983. 

[2] I. Antoniadis, C Kounnas and R. Lacaze, Nucl. Phys. B211:216, 1983. 

[3] G. Rossi, Phys. Rev. D29:852, 1984. 

[4] C. Coriano' and C. Savkh, Comput. Phys. Comm. 118:236,1999. 

[5] W. Furmanski and R. Petronzio, Nucl. Phys. B195:237, 1982. 

[6] M. Botje, QCDNUM16: A fast QCD evolution program, ZEUSS Note 97-066. 



24 

[7] A. Chuvakin and J. Smith, |hep-ph/991150l . 

[8] J.H. Da Luz Vieira and J.K. Storrow Z.Phys.C51:241-258,1991. 

[9] L.E.Gordon and G.P. Ramsey Phys.Rev.D59:074018,1999. 

[10] H.L. Lai et al, Phys. Rev. D55:1280, 1997. Phys. Rev. D51:4763, 1996. 

[11] J. Botts and J. Blumlein Phys.Lett.B325:190-196,1994. 

[12] L. Clavelli, |hep-ph/9812340| . 



25 



4.0 



3.0 i 



2.0 \ 



\ 



up 

down 

strange 

gluon 




Figure 1: u, d, s and gluon distributions at Qi = 5 GeV 
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Figure 2: xu{x) and xd{x) evaluated with Qi = 2.0 GeV and Qf = 100 GeV with mx = 10 
GeV in the standard (reg) and susy evolution 
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Figure 3: Gluon distributions with Qi = 2.0 GeV and Qf 
GeV. The regular and the susy evolution are shown 
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Figure 4: xq^~^\x) (singlet) quark distribution evaluated with Qi = 5.0 GeV and Qf 
GeV with m\ = 10 GeV in the standard (non-susy) and susy evolution 
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Figure 5: Comparison of the supersymmetric gluon and gluino distributions for mx = 30 GeV 
and Qf = 100 GeV. 
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Figure 6: Gluon distribution for 3 values of the gluino mass 100, 150 and 200 GeV and a final 
evolution scale Qf = 1 TeV. Shown is also the regular (reg) evolution of the same distribution 
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Figure 7: Gluino distribution for Qf = 1,2 and 5 TeV and rrix = 100 GeV 
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8: Gluino and gluon distributions for very large final evolution scales Qf = 5 and 10 
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Figure 9: Gluino distributions for a varying rrix (40, 100, 200 and 250 GeV) with a fixed final 
evolution scale Qf = 1 TeV. 
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Figure 10: u, d quark distributions rrix = 250 GeV with a varying final evolution scale Qf = 5 
and 10 TeV. Shown is also the strange quark distribution with Qf = 5 TeV. 
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Figure 11: Dependence of the total 2-gluino cross section cTpp^w on the factorization scale in 
the QCD and SQCD cases {mx = 250 GeV). Shown are the factorization scales Qfact= 510, 
600 and 700 GeV. We show both the QCD and the SQCD resuhs 



32 







II 


' 1 ' 1 ' 1 ' !__'_, 






QCD(510GeV) 




10-^ 


— 


■■■■ QCD(700GeV) 






~ 


-- SQCD(510GeV) 


^"^ .^'^ - 




- 


■- SQCD(700GeV) 


^ ^ 
^ 




^ y H 




^^ / 




• . 




/' / 




X^ ./ 




/ / 


'>? -4 


/ ./ 


^ 10 


/ / 


X 


/ 




/ / : 




/ /' ^ " 




'' / ^..^'^ 








// ^y^ 




/ / y^ 


_s 


1 ' y^ 


10 






- ' ! y\' - 




- 1 ' /\' - 




1 /I.I 1 1 1 \/ .• 1 1 1 1 1 1 



600 



800 1000 1200 1400 1600 1800 2000 

X 



Figure 12: Dependence of the total 2-gluino cross section Ogg^w on the factorization scale in 
the QCD and SQCD cases {m\ = 250 GeV). Shown are the factorization scales Qfact= 510 and 
700 GeV.We show both the QCD and the SQCD results. 
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Figure 13: Dependence of the total 2-gluino cross section o-^g-^AA on the factorization scale in 
the QCD and SQCD cases {m\ = 250 GeV). Shown are the factorization scales Qfact= 510 and 
700 GeV.We show both the QCD and the SQCD results. 
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1 Introduction 



The study of the QCD scahng violations is an important chapter of high energy physics and a 
very important tool for the analysis of future data at the LHC At large energy and momentum 
transfers, the underlying quark-gluon dynamics is light-cone dominated and controlled by a 
mechanism of collinear radiative emission described by logarithmic corrections to the lowest 
order cross section. In this picture, based on the parton model, initial quark states are assumed 
to be massless and the running of the coupling is linked to the number of flavours nj included 
in the evolution. Crossing intermediate thresholds opens up new channels and new dynamics. 
There is a widely used formalism of perturbative QCD that we think is worth to extend to the 
super symmetric case and which might be useful for experimental searches of supersymmetry at 
the LHC. The picture is particularly appealing if a "supersymmetric content " of the proton 
is found. However, predictions tied to this picture may be used to rule out a possible light 
gluino and/or a light squark. Other applications involve the study of nucleon collisions with 
the atmosphere as in Ultra High Energy Cosmic Rays, where the center of mass energy of the 
primaries can reach several hundreds of TeV's. 

This formalism involves the study of scaling violations induced on the initial state by the 
opening of supersymmetric channels prior to reaching the hard scattering phase. The analysis 
requires the notion of supersymmetric parton distributions, supersymmetric factorization for- 
mulas for supersymmetric DIS and (susy) hadron-hadron collisions on which we elaborate in 
detail. We remark that supersymmetric versions of parton distributions are easy to define, by 
a natural generalization of the usual DIS approach. The phenomenological validity of these 
extensions are clearly linked to the possibility of detecting scaling violations in precision mea- 
surements of final states in p-p collisions, for instance in Drell Yan, or in other processes with 
a distinct final state. With this analysis in perspective we start developing accurate tools that 
we will use in a separate work for a quantification of rates for gluino/squarks effects at the 
LHC. Although there is no evidence of supersymmetry at current energies, or of supersymme- 
try coming from the initial state, quantifying with accuracy these effects we believe, is still an 
interesting task. These studies are also of direct theoretical relevance, since they tell us more 
specifically how to merge the usual parton model dynamics with the new elementary states 
that supersymmetry predicts. 

In this paper we start analizing these issues and make a first step toward quantifying the 
impact of supersymmetry in the initial state. Our analysis, in this paper, is focused on the 
supersymmetric DGLAP evolution, which is more involved compared to ordinary QCD. We 
solve the equations using a recursive algorithm that we formulate and that we have tested 
which allows to perform a direct match between the various regions of the evolution as we 
increase the factorization scale and allow for new (supersymmetric) channels. Part of the 
analysis is of technical nature and several sections deal with the implementation of the method. 
We then study in 2 final sections the implications of the susy evolution and compare it to the 
standard QCD one. Other aspects of the evolution, including some applications for collider 
processes will be analized in a companion paper. 



2 The Method 

Before moving to Supersymmetric QCD (SQCD), we briefly illustrate the method as it applies 
to the case of ordinary QCD. This simpler case will help us establish notations that will be 
later extended to the supersymmetric evolution. We build on previous work of Rossi [3] , which 
has been used before in the previous literature [9, 8], although never fully documented in its 
code implementation. We combine that method, originally suggested for QCD, with recent 
approaches that use an analytic evaluation of the convolution integrals ("weights") [6], and 
generalize it to SQCD. A recently application to QCD in next-to-leading order (NLO) of this 
approach has been implemented by Chuvakin and Smith [7] in their analysis of an evolution 
scheme with a varying number of flavours. To make the discussion self-contained, we have 
elaborated in some detail on the method in an appendix. 

The two-loop running of the coupling constant is defined by 

«(Q^)_2 1 f^ A lnln(gVA^) ^ ^^ 1 ^^^ 



27r /3oln(gVA2)V (3o HQV^^) Mn^(Q2/A2) 

where 



Po = yCc - -TRHf 



A = f C^ - fCaUf - 2Cpnf, (2) 



where 



iV^ _ 1 1 

Cg = N, Cp = ^^^. Tr=- (3) 

and N is the number of colours while nj is the number of flavours. 
The solution for the running coupling is given by 

ait) = ^e-^°/2* (4) 

with <y{Ql) = q;(0), and Qo denoting the initial scale at which the evolution starts. The 
evolution equations are of the form 



.2 d , „2N "(Q^^ 



^^^^^'^'^^ " ^^P(-)(x,a(g^))®X.(a:,g'), 



(5) 



with 



X^{x,Q') = q^^'-\x,Q') - — gW(^,Q') (6) 



Up 



for the non-singlet distributions and 






i?(±)(x) = PS(:r)-£^Pf(x) 



(7) 



for the singlet sector. 

We have defined, as usual 

rif 

Qi'^ = <li- Qi, Qi^^ =qi + Qi, g^^^ = ^ = Y1 (li^^- (8) 

i=l 

We introduce the evolution variable 

Po a(<5o) 
which replaces Q^. The evolution equations are then rewritten in the form 

j^qt\t,x) = (^P(0)(x) + ^i?(_)(a:) + ...)®gj-)(t,a:) (10) 

Q'J^X^{X,Q') = (p^'^\x) + ^R^^)ix)\0X^ix,Q'), (11) 

In the new variable t, the kernels of the evolution take the form 



Rix) = P«(x)-Ap(o)(^). (13) 

Equations (10) and (11) are solved independently for the variables g| and Xi respectively. 
Finally, the solution q^~^^ of eq. (12) (or the singlet equation) is substitued into Xi ^^ order to 
obtain q} ' . 

The equations can be written down in terms of two singlet evolution operators E±{t, x) and 
initial conditions q±{x,t = 0) = q±{x) as 

j^E^ = P±^E±, (14) 

whose solutions are given by 

qt\t,x) = E^_)0qt^ 

Xi{t,x) = E(^+)®Xiix). (15) 

The singlet evolution for the matrix operator E{x, t) 



EGF{x,t) EGG{x,t) 



(16) 



'-^ = P^Ei.,t) (17) 

is solved similarly as 



G{t,x) -Ht,x}(S\ ^^^^ 



The unpolarized leading order kernels are expanded in a as 



P^^i^,C.s) = (^) PS\x) + (^)'4'^(^) + - (19) 



Their LO expressions are given by 

^'"--^■.(ii^)^ (20) 

for the non-singlet sector, and by 



p(0)/ \ _ p{0) 
^qq K-^) — ^qq,NS 



PjJ)(x)=2T^n;(x2 + (l-x 

p(o)(.)=C.i±il^^ 

P(;)(x) = 2iV, ((Y3^ + \-2 + x{l-x)^+ ^5{l - x) (21) 

in the singlet sector. Some simple identities for the "plus' (+) distributions and their definition 
are given in the appendix. 

In the actual numerical solution of the equation, one would like to have at hand a recursion 
relation which can be implemented in a computer program dynamically at run-time. In order 
to develope recursion relations we proceed as follow. We first observe that as far as we are not 
resumming double logarithms in Q and x, the solution eq. (18) can be expanded as a series of 
convolution products 



q{x, t) = qo{x) + tP ® qo{x) + -P ® P ^ qo{x) + ... (22) 



and transformed into a recursion relation for some coefficients Cn{x) 



Co{x) 
Cn+1 



6{x-l) 
P ® Cn{x) 



equivalent to the expansion for the evolution operator 

n! 



E{x,t) = J2CNix)- 

n=0 



(23) 



(24) 



Written as a recursion relation, the equation is as easy to implement as a calculation in moment 
space, but with no need to perform any moment inversion. The possibility of using this expan- 
sion as an alternative way to evolve parton distributions fast and with accuracy especially for 
larger sets of coupled equations, such as for supersymmetric theories, is the subject of the follow- 
ing sections. Specifically, Rossi's ansatz [3], originally introduced in the context of the photon 
structure funtion, differs from (24) only by an over all normalization Cn = (~1)^ (2//3o) 
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(25) 



Connecting this expansion to other expansions is also pretty straightforward. A very elegant 
basis in which to expand is the Laguerre basis, introduced by Furmanski and Petronzio [5]. A 
relation between Rossi's basis and the FP basis can also be derived but the result is not 
particularly illuminating. 

Inserting the expansion (25) into the evolution equation one gets some recursion relations 
for the functions An+i{x) in terms of the An{x). These are obtained by comparing left hand 
side and right hand side of the evolution equations after equating the logarithmic powers with 
a running strong coupling constant a{Q'^). A pretty detailed study in the polarized case can 
be found in ref. [9]. For computational purposes, the recursion relations for the evaluation of 
the An{x)^s can be implemented in various ways. The one that we have implemented involves 
the direct solution of the recursion relations as illustrated below. In leading order in q;(Q^) in 
QCD we get for the unpolarized kernels 
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int the non-singlet case and 
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Po Jx y [z 

in the singlet. Extensions to the NLO case of this approach are pretty straightforward, and the 
procedure will be more clear once we will discuss its implementation to solve the Supersym- 
metric DGLAP (SDGLAP) equations. 

3 Supersymmetric scaling violations 

In A^ = 1 QCD gluons have partners called gluinos (here denoted by A) and left- and right- 
handed quarks have complex scalar partners (squarks) which we denote as ql and qn with 
q = Ql + Qr (for left-handed and right-handed squarks respectively). 

The interaction between the elementary fields are described by the SU{3) color gauge in- 
variant and supersymmetric lagrangean 



+qj Pq, + D^qnD'^q'' + D^q^D'^q'^ + tgV2 (x-RqlT'^qu + ^'^LqlnT'^qm - h. c 



■^9^[qL^T\Li-q\i^T''qli,j +mass terms 



(28) 



where a runs over the adjoint of the color group and i denotes the number of flavours over 
which we sum. 

Past studies of these effects lead to the conclusion that the information that supersymmetric 
initial states carry along the evolution can be easily absorbed into scaling violations coming from 
ordinary QCD. These studies require the knowledge of the matrix of the anomalous dimensions 
(for all the moments), or of the corresponding Altarelli-Parisi (DGLAP) kernels. We recall that 
the leading and next-to-leading anomalous dimensions are known [2] since long ago, at least in 
the case of partial supersymmetry breaking, in a scenario characterized by a light gluino and 
a decoupled scalar quark (we present in an appendix the leading order form of these kernels 
for completeness). In this work we consider the possibility of a radiatively generated gluino 
distributions, in analogy with the case of standard QCD for the gluons. In fact, in QCD one 
can use a simple model of the distributions at low Q {Q = Qo) and run the DGLAP equations 
up in energy in order to generate distributions of gluons at any higher scale Qf. Of course, it 
is well known that in QCD such initial evolution scale is model dependent. Various amplitudes 
at higher energy emerge, depending on the underlying assumptions on the form of their initial 
shapes. 



This approach has been widely used in the hterature and can be a way to connect low 
energy quark models -which have no gluons but a phenomenological confining potential - to 
true QCD scattering amplitudes. One can assume a zero distribution of glue at the lowest scale, 
or a model dependent non zero one, and then use the evolution equations to dress the matrix 
elements describing the distributions by logarithmic corrections. Once the final scale -here 
identified as the factorization scale of the process - is reached, the distributions are convoluted 
with the usual -on shell- parton cross sections to generate the full hadronic cross section. 

In hadronic collisions, Qf is usually a fraction of the center-of-mass energy, or a fraction 
of the large pt of the final state jets. It is not uniquely identified and a different choice of Qf 
underscores a drastic sensitivity of the expansion on this scale. 

Alternative choices for density of the constituents at low Q are also possible, but, ultimately, 
whatever the choice for the parton structure, it has to match the scattering data available from 
DIS and colliders experiments. In the case of exact SQCD, it is natural to parametrize the 
parton distributions, now with a susy content in the form 

F2{x,Q^) = xe^i{qi{x,Q^)+qi{x,Q^) + qiii{x,Q^)+qiji{x,Q^) 

q^L{x,Q') + lL{x,Q')), (29) 



where the new elementary constituents (squarks and gluons) can be vanishing at low Q and be 
radiatively generated by the evolution, similarly to the gluon case (in standard QCD) contem- 
plated above. In the case of a broken susy, the form of F2 does not change from the standard 
QCD form. 

A gluino or a squark parton distribution (a more detailed analysis will be presented else- 
where) is the exact correspondent of a gluon or a quark parton distribution, namely a light- 
cone dominated diagonal (non local) correlation function in spacetime, Fourier transformed to 
(Bjorken) x-space. 

Similarly to the QCD case, in the case of exact A^ = 1 supersymmetry we define singlet and 
non-singlet distributions 



g^(x,Q^) = ^(gi(x,Q2) -gi(a:,Q2) 
1=1 

qv{x,Q^) = ^(gi(a:,Q2) -gi(x,Q2) 

i=l 

q+{x,Q') = ^(g,(x,Q2) + g,(x,Q^) 
1=1 



1=1 



The evolution equations can be separated in two non-singlet sectors and a singlet one. The 
non-singlet are 
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and the singlet, which mix qy and qy with the gluons and the gluinos are 
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(32) 



There are simple ways to calculate the kernel of the SDGLAP evolution by a simple extension 
of the usual methods. The changes are primarily due to color factors. There are also some 
basic supersymmetric relations which have to be satisfied that will be analized below. They 
are generally broken in the case of decoupling. We recall that the supersymmetric version of 
the P function is given at two-loop level by 
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(33) 



with Cx = Ca = Nc for the case of Majorana gluinos and the ordinary running of the coupling 
is replaced by its supersymmetric running 
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(34) 



The kernels are modified both in their coupling {a -^ a'^) and in their internal structure 
(Casimirs, color factors, etc.) when moving from the QCD case to the SQCD case. In our 
conventions an index "S" stands for a supersymmetric component (regular, i.e. non supersym- 
metric, kernels do not carry such an index), but we will omit it when obvious. 



4 Evolution and Matching 

Susy is necessarily broken in the real world and therefore, the way the breaking occurs dictates 
both the mass spectrum and guides the structure of the QCD scaling violations as well. There 
are several parameters that appear in the evolution, mx and m^, the masses of the gluino and 
squarks, hence a complete analysis includes various scenarios, on which we briefly elaborate. 
In a realistic scenario with a broken susy, the squarks have much larger mass compared to the 
quarks and the gluinos have a Dirac or Majorana mass mx (or a combination them). In our 
case, in order to establish the evolution scales which are of phenomenological interest, it is 
convenient to assume 1) that the scalar quark decouples from the evolution and 2) that the 
2-gluino production threshold m2x = 2mx is the intermediate scale, separating in Q^ the regular 
QCD region from the supersymmetric one. In this scenario some of the splitting functions - 
specifically Pgx and Pxq - are zero, since no collinear emission is associated with massive partons, 
unless the symmetry is effectively restored. 



We impose the separation condition 

Qo < "^2A < Q 
and only one non singlet evolution equation is considered 
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The simphfied singlet equations, which mix qy and qy with the gluons and the gluinos 
become 
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Formally, we can write down the expression of the complete solution in the form 
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for the non-singlet equation and 
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for the singlet one. In analogy with eq. (25), now we introduce supersymmetric coefficients 
A^i^x) beside the usual non supersymmetric An{x), and impose recursion relations on the initial 
ansatz q{x, Ql) of the form 
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n = 1, 2, ..., which solve the equation in the non-singlet sector. 
In the singlet sector we get the solution 



G{x, Qq) 
\{x,Ql) 



G{x, Qq) 


qix,Ql) 



E 

n=l 



n! 






+ E 



ra=l 



1 



AUx,Q'] 



log" 



a(m2x] 



constructed recursively through the relations 
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(4J 



n = 1, 2, ..., having used a vanishing gluon density at the starting point of the supersymmetric 
evolution. 



5 Supersymmetric Relations 

An exact supersymmetric scenario is probably interesting only for analizing theoretical issues 
concerning the evolution, or to study the shape of the distributions at extremely high energies, 
when, again, all the supersymmetric partners effectively become mass degenerate. This scenario 
can take place at few TeV's or at several TeV's, depending upon the assumptions underlying 
the way supersymmetry is restored. Here we simply focus our analysis on a light, up to an 
intermediate-mass Majorana gluino. Other aspects of the evolution, such as the interplay 
of Dirac and Majorana gluinos and the impact of various patterns of susy breaking will be 
considered elsewhere. 

As we have already mentioned, the study of the anomalous dimensions - in leading order - 
for A^ = 1 QCD can be found in [1] and now we are going to elaborate on that. 



We recall that in a regularization scheme which is manifestly supersymmetric there are some 
super symmetric relations which are satisfied by the anomalous dimensions. This is true, for 
instance, in the Dimensional Reduction DR scheme. We recall that to leading order the MS 
scheme and the supersymmetric DR scheme give coincident results. We recall that in the DR 
scheme the traces are kept in 4 dimensions, and the chiral projectors are treated as usual, with 
a completely anticommuting 75. The loop momenta are evaluated in n-dimensions. In other 
schemes, such as the t'Hooft-Veltman scheme, 75 is instead partially anticommuting. The MS 
scheme for chiral states is usually based on this second definition. The relation between the 
two schemes is to leading order 

p(0) _ p{0) .^gX 

MS DR ^^^ 

valid both for the non singlet and the singlet anomalous dimensions. 

The same is true also for the factorization scheme dependence of the coefficient functions. 
As for the coupling constant, the definition of a in the two schemes, the expression of aj^g can 
be used also in the DR scheme as far as the Aqcd scales of the two schemes are related by 
A-^^ = Aj^exp (Ca/QPo)- This last change is tiny and will be neglected. 

Part of the supersymmetric kernels can be obtained at this order by a simple change of 
group factors, from the standard QCD kernels. For instance, the substitutions Cp —^ C\ 
{Cp = {N^ — l)/{2Nc)) and nj -^ nx, allow us to obtain the kernels {\\, \g,gX) from the 
ordinary kernels {qq,qg,gq)- We use Cx = Ca = N^, the number of colors. As for the choice 
of the type of representation for the gluino ( Majorana or Dirac), we recall that in the Dirac 
case we set nx = 2Ca while in the Majorana case nx = Ca- The organization of the terms in 
the splitting functions may differ from reference to reference, due to the various manipulations 
one can perform on the plus distributions. 

The supersymmetry relations are given by 

(50) 
and by 
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(51) 

In general, for a decoupled scalar quark, the only symmetry one would expect is eq. (50). In 
the case of Majorana gluinos, for a scenario with a decoupled squark, it is interesting to observe 
that this relation remains valid for a; < 1 as well. It can be extrapolated to include the x = 1 
point in the case of zero number of flavours (supersymmetric gluondynamics) . 

The evolution has to respect - both in the case of exact susy and of susy breaking - 1) baryon 
number conservation and 2) momentum conservation. There are two sum rules associated 
with these conserved quantities, which are generally used to constrain the phenomenological 
parametrizations of the distributions in direct applications. Below the ni2x threshold we need 
to satisfy the usual QCD relations 

dx {xG{x, Q^) + a;g(+) (x, Q^)) = 1 (52) 



for momentum conservation and 
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Jo 

for baryon number conservation. These require that the anomalous dimensions satisfy the 

relations 

1 
xdx {Pgi{x) + Pgi{x)) = 

/ dxPNsix) = 0, 
Jo 

(54) 
respectively, where i = q,g. In leading order the second equation is simply 

f dxPf^{x) = 0. (55) 

Moving above the 2-gluino threshold the momentum sum rule becomes 
"1 



dx {xG{x) + x\{x) + xq^^'{x)\ = 1 (56) 

for momentum conservation, while eq. (53) remains unaltered. We get the new momentum sum 
rule (or second moment sum rule) 

f xdx (Pgiix) + Pg,{x) + Pxi) = 0, (57) 

with i = q,g, A. The supersymmetric version of eq. (55) is simply obtained by replacing P^g\x) 
by it supersymmetric counterpart P^J^^^x). It can be checked easily that the kernels in the 
appendix satisy these relations. 

In the case of exact supersymmetry we need to keep into account the scalar quark contri- 
bution in both equations. In particular, the equation for the second moment becomes 

/■I 

/ xdx (xG{x) + x\{x) + xq^^'{x) + xq^'^\x)) = 1, (58) 

which implies that 

fxdx (Pgiix) + Pgi{x) + Pxi + Pgi) = 0, (59) 

with i = q,g, A, q. As for baryon number conservation, equation (53) gets modified into 

/ dx {q^-\x) + q^-\x)) = 3 (60) 

and using (31) one gets two relations 

^' rfa: (Pf, + P/,) = 



(61) 



(below we will drop the susy index S in front of the kernels when obvious). The first relation 
in the equation above, for instance, clearly implies that the end point contributions in the 
ordinary Pgg kernel are to be modified in order to insure conservation of baryon number. In 
the case of a susy breaking scenario such additional end-point contribution is trivially absent 
and the form of the Pgg kernel remains the same, except for the replacement of a with o;"^, the 
super symmetric coupling and Pq with Pq . 

6 Applications 

In order to apply the formalism to SQCD we proceed by discussing the method for a scenario 
with a broken susy (decoupled squark). We proceed from the non singlet case, analizing in 
details the intermediate steps of the algorithm. 

We start with a solution of the standard non-singlet DGLAP equation assuming the bound- 
ary condition Aq (x) = q^{x, Ql) and run the equation up to the gluino threshold and construct 
the coefficients recursively. We arrest the coefficient up to a desired order (n), which can be as 
large as 30. In general, the rate of convergence of the asymptotic expansion changes with the 
value of the momentum Q. 

The solution is then constructed in the region Qo < Q < m2\ as 

/(..-.) = E^luf^) (62) 

and used as initial condition for the next stage of the evolution, which involves the region 
'^2A < Q < Qf, with Qf being the final evolution scale. 

At the next stage, we set Aq{x) = 6{1 — a;) q^{x,m2\) and solve recursively using the 
super symmetric version of the kernels. The strong coupling constant Q;((5^) and its running 
are replaced by their supersymmetric version a'^(Q^), and so are the coefficients of the beta 
function {(3i -^ f3?). Finally the solution is written down in the form 



x: ^^ log" f^^^!^] + 1^ ^ log" ( ^i^^>^A 



n=0 



q'{x,Q') = E ^^ log" RTTf^ +E^^ log" RtT^ (63) 



n! V "('^o) J ~i n\ \ a{Q) J 



where n' is the index at which we arrest the supersymmetric recursion. In our implementa- 
tion we have kept the values of n and n' very close. 

In the singlet case the procedure is not much different. We evolve according to the standard 
DGLAP equation in the region below the 2-gluino threshold, having set to zero the contribution 
from the gluino at the beginning (up to the supersymmetric threshold). This is equivalent to 
having set A^{x) = for any n, which means that in the region below m2\ there is no radiative 
production of gluinos in the initial stage. We iterate the recursion relations up to a given value 
n of the index n with the standard DGLAP. The initial conditions for the next stage of the 
evolution are then fixed by the relations 

Al'^{x) = q^{x,m2x), 

= i:^>og"(%4). (64) 

Af^{x) = G{x,m2x) 



= t^los-i^] (65) 

(66) 
After this, we determine recursively the coefficients An of the supersymmetric expansion 

""^^ P^ Jx y y-x 

Po Po "'^ 1/ Po 

Pq Jx y ^ 

(67) 



4A .X _ ,C'a /■irft/ |/A^(|/)-a:A^(a:) 



2§r^(l + ^Mn(^)-4c',A^.(^) 
Po "'^ 1/ Po 

4nA /' ^(1 - 2^ + 2z^))i^(y) - 4§i^(x) log(l - x), (6J 

Po "'^ ^ Po 



AU. = -^cj'^(--2 + z)AUy)-^cj'^r--2 + z)Atiy) 



^^ ^ ^ y^nKyj-^^^nK-^J _ 4^ i^g(i _ ^)^, (^) 



/^o* Jx y \z J " /^o* Jx y V2; 

, CU /•! rft/ yi^(t/)-a:i^(x) _ ^C;^ 
' Po Jx y y-x Pi 

Pq Jx y \z / 



Finally, we construct the solution in the form 

^0 n! V "('^o) / ^1 n! \a[m2x)) 



(69) 



(70) 



(71) 



where we have arrested the supersymmetric recursion up to the index n' . Here f{x,Q'^) indi- 
cates a singlet quark, a gluon, or a gluino distribution, with A{(x) and A{(x) denoting their 
corresponding coefficients in the expansion. Singularities emerging from the lower integration 
point {y = x) are the tricky part of the game, as expected, but can be handled with various 
techniques. They will be discussed briefly in the last section. 



7 NLO Extensions 

Let's now move to a next-to-leading order (NLO) analysis of the evolution. Our discussion here 
is partial and does not include contributions due to the emergence of new anomalous dimensions 
as we move across the supersymmetric threshold. In recent work [7] it has been shown that in 
the ordinary distributions of quarks and gluons of QCD these effects are important, especially 
at small-x. These changes are expected to produce only a slight modification of the algorithm 
presented below, and simply amount to a modification of the boundary condition as we move 
across the 7n2\ point. They will not be analized any further in this work and will be assumed 
to be negligible. As a second point, we remark that the extension of the procedure outlined 
below is easy to generalize to the more general case of exact susy. 
To NLO the ansatz becomes 

and inserting the usual running of the coupling 

p{a) = — -;— tt — „ q: (73) 



c^log(Q2) ^^ ^ 47r 167r2 

we get the recursion relations 



2 

An+l = - — An{x) 



A 







Bn+i[x) = -Bjx 



'r,{x) - (^A^^.ix)] - -^P(^) ® A„(x) - ^P^'') ® B^x) 



4/3o ""^^ V AttPo "'' Po 



V^Po / 47r/;/o Po 

(74) 

which are solved with the initial condition Bq{x) = 0. The initial condition for the An{x) 
coefficients (i.e. Aq{x)), is specified as in the previous section, with q{x,Ql) identified as the 
leading order ansatz for the initial distribution, i.e. 

Ao{x) = 6{1 - x) ® q{x, Ql) = g(^°)(x, Ql) (75) 

Running the RGE's below the gluino threshold region (up to Q = m2\) and solving the 
recursion relations (74), we get the solution (arrested at a recursive index n) 

which is used to fix the intial condition for the second stage of the evolution, the supersymmetric 



one 



q{x,m2x) = q {x , m2x) + a{m2x) q {x,m2x), (77) 



with 

The supersymmetric recursion relations are then given by 

Ao{x) = 6{l-x)(g)q^^{x,m2x), 
Bo{x) = g^^^(a;,m2A), 



Po 



A ;; /„\ \ ^ r.a)5 ^ 4 /„\ 2 



5„+i(x) = -5„(x)-(^^i„+i(x)j-^P«^®i„(x)-^pW ®P„(x) 



^JS' 






(x)) - -i^p(i)^ ® i„(x) - -|p(o)5 ^ ^^(^)_ 
/ ^Po Po 

(79) 



We finally construct the general solution in the form 
Eq. (79) can be expanded in components, since it is valid in matrix form 

B'M^) = -^r(:^)+-C(n^r®^n^(^)+n^r®^n(^)+n^?'®^n(^) 

2Po 

Po 
and similarly for the evolution of the gluon density 

2Po 
4 Wr ® ^f (^) + P^'J'® ^Ux) + P«" ® i^(x 



4/35' 

Po 



(82) 



while the gluino density is obtained using the recursion relations 



Po 



(83) 



Notice that as initial condition for the gluino distributions we take an identically vanishing 
function at the scale Q = m2\ both in leading and in next-to-leading order. 

Ao{x) = 
Boix) = 0. 

(84) 

8 Numerical Results 

It is expected that a large gluino mass, for a fixed factorization scale Qf in the evolution, 
lowers the size of the scaling violations and their impact on the supersymmetric cross section. 
On the other side, scaling violations induced by the susy evolution should grow as we raise 
the final evolution scale. Therefore it seems natural to study the effect of the susy evolution 
in two different setups 1) for fixed mx and a varying Qf or 2) for a varying mx at a given 
factorization scale Qf. We have performed both studies and the results are shown in figs 1-14. 
The implementation of the unpolarized first stage (QCD) evolution is performed in the MS 
scheme, which is by now standard in most of the high energy physics applications. We introduce 
valence quark distributions qvi^^Qo) and gluon distributions G{x,Ql) at the input scale Qo, 
taken from the CTEQ3M parametrization [10] 

-x)'^^{l + A3x'^^). (85) 



q{x) 


= AqX 


Specifically 






xuv{x) 




xdv{x) 




xG{x) 


.^(^^ 


1 + d{x) 



1.37x°-^^^(l-x)3-7^[l + 6.25x°-H 
0.801x°-^9^(l-x)^-i9[l + 1.69x°-37^] 
0.738x-°-=^^^(l - xf-^^[l + 7.30a;] 

0.547x-°-286(l _ a;)8-34[l + 17.53.] 

, , u(x) + d(x) 

xs{x) = 0.53; 2 

x{d-u) = 0.75x^-^^(1 -a;)^-3^(l + 30.0 a:) (86) 

and a vanishing anti-strange contribution at the input. In figs. 1-13, where we have studied the 
valence quarks, the gluon and the gluino distribution for a varying gluino mass ranging from a 



light to an intermediate gluino (10-40 GeV) up to a value of 250 GeV. These results generally 
point toward small scaling violations, which become more appreciable as we move closer the 
smaller-x region (in particular for gluons and gluinos). We have chosen as initial evolving scale 
Qi = 2 GeV in all the runs. 

Fig. 1 shows the shape of the distributions close to their initial evolution scale, with Qi = 2 
GeV and Qf = 5 GeV. In fig. 2 we show the scaling violations induced by the susy evolution 
on the u and d quark distributions, from the initial scale up to the final scale of 100 GeV for a 
gluino mass of 10 GeV. In this figure we show the regular versus the supersymmetric evolution 
for these two distributions. The modifications appear to be quite sizeable. The presence of 
supersymmetry in the evolution shows up as a lowering of the maxima of the distributions with 
a shift toward the small-x region. 

A better distinction between the non susy from the susy result is illustrated in fig. 3, 
which shows a comparison between regular and susy evolution for gluons. As expected, these 
differences get more pronounced moving toward the region of smaller x, due to the rise of the 
gluon distribution and to the small-x structure of the kernels. In fig. 4 we show the singlet 
q^^^ distribution. The difference between the regular and the supersymmetric evolution is slim, 
given the low factorization scale (100 GeV) chosen in this run. 

Gluino and gluon distributions are shown in fig. 5. The gluino distribution is approximately 
2 orders of magnitude smaller compared to the the supersymmetric gluon distribution and 
grows fast at small-x. Fig. 6 shows a plot of the gluon distribution for a large final factorization 
scale Qf = 1 TeV and a varying gluino mass (100,150 and 200 GeV respectively). Shown is 
also the regular evolution of the gluon density, which is lower than the susy ones at larger x 
vales and faster growing at smaller x values. Fig. 7 shows the factorization scale dependence 
of the gluino distribution for a sizeable gluino mass (100 GeV) and large factorization scales 
Qf = 1,2 and 5 TeV. We plot in fig. 8 the gluon and the gluino distributions for very large 
evolution scales for a realistic gluino mass of 250 GeV and varying factorization scales. We 
have chosen Qf = 5 and 10 TeV respectively. The dependence on the final scale appears to 
be quite mild. Fig. 9 shows the variation of the gluino distribution for different gluino masses 
(40,100,200 and 250 GeV) and a fixed final evolution scale Qf = 1 TeV. There is a reduction of 
the small-x growth of this distribution at smaller-x values as the mass of the gluino is raised. In 
fig. 10 we show the dependence of all the quark distributions on the factorization scale Qf = 5 
and 10 TeV in the case of a supersymmetric scale m2\ of 250 GeV. 

We study the impact of these corrections on future collider experiments by showing results 
for the 2-gluino production in a p-p collision. As can be seen from fig. 11, the production cross 
section is quite small -compared to standard QCD rates-, but gets enhanced by the inclusion 
of susy scaling violations especially at larger energies. We have set the gluino mass at 250 
GeV. In the same figure we show the dependence of the cross section on the factorization scale 
(510,600 and 700 GeV). The dependence is sizeable, although the total rates remain small 
compared to the QCD background at these energies ( 650 GeV - 2 TeV). We have shown in 
figs. 12 and 13 the partial contributions to fig. 11 coming from the g g ^ XX channel and the 
qq channel. The gluon contribution is the dominant part in both the regular QCD case and in 
the supersymmetric case. 



9 Conclusions 

We have solved the supersymmetric DGLAP equation in a scenario with a coupled gluino and 
a decoupled squark using an algorithm based on x-space. In particular, we have illustrated the 
evolution of the distributions of quarks and gluons and their supersymmetric versions using a 
radiative model. Although the window on a light gluino is now rapidly closing and the hope 
to detect supersymmetry from scaling violations in the initial state at the LHC with a light 
gluino is slim, the possibility of analizing experimentally the impact of heavier supersymmetric 
particles in the initial state remains, however, an important issue. In a specific example (p- 
p — > AA), we have seen that the total cross section is sensitive to initial state susy scaling 
violations and on the factorization scale chosen. For 2-gluino production, for instance, the rate 
is small, much below the usual QCD background, but gets sizeably enhanced when susy is 
included. We believe that it is interesting to see how much we can rise the gluino mass and still 
obtain a signal on a final state which can not be compensated by the usual re-adjustement of 
the several parameters that describe the ordinary QCD parton distributions, or is comparable 
to it. In particular, the strong factorization scale dependence may be reduced by the inclusion 
of radiative correctisons in the initial state. This and other related issues will be analized 
elsewhere. 
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9.1 Appendix 1: The Kernels 

Various relations among different types of "+" distributions can be derived. The kernels given 
in the literature may differ in their final expressions due to rearrangements of the correspond- 
ing "+" functions. We give here the expression of these kernels and illustrate some of the 
manipulations needed to reorganize them in a standard form. 

In the derivation of the recursion relations we have used the identity 

'^^ ' V„(i,)^/'^ ''-'-''''-"^"'"' +log(l-.M„M (87) 



Ix y \il -x/y)+) Jx y y-x 

Similarly, it is not hard to show the identity of the convolution products 



'^An{x) = / — —An{y)-A{x) dz- A{x) 



xj Jx y{y-x)y Jo 1 - z Jx y y^{y - x 



with 



M-l-x + |*(l-x) 



A = C ^y^ ^^""^^^ -xAn{x) 



y y-x 



3 /■! dv 

- 2log{l-x)A^ix) + -Aix)- -Ji{l + z)ArXy) 

2 Jx y 

(89) 



After some manipulations one can show that the two expressions given above are equal and 
therefore 



'l+x^\ 2 ^ 3^,^ 

=- ^ 1 — X -\ — oil — x) 

l-x j . {l-x)+ 2 ^ ' 



(90) 



In leading order, the supersymmetric expression of the standard (qq,qg,gq,gg) QCD kernels 
are obtained by replacing a -^ a"^, and the dependence on the coefficient of the j3 function by 
its supersymmetric version 13q — > /3^. 

The leading order expression of the kernels in the case of a decoupled scalar quark are given 
below. We have omitted the superscrit "S" for simplicity. We remark that the usual QCD 
kernels, after the embedding in the supersymmetric evolution, do not acquire new terms at 
a; = 1, except for Pgg. 



p(0) 
99 



>(0) 



SA 



P(0) 

9q 

p(0) 
p(0) 



p(0) 

p(0) 

p(0) 
19 

p(0) 

p(0)5 
11 



2Ca 



1 



1 



(1 — x)^ X 
l + (l-x)^' 



x(l — x] 



f^a-a:) 



X 



Pi°) = C, 



1 - 2a; + 2a;' 

2 

(1-x)^ ~ 



99 

nx 
Cx 




^? = % 
nf{l-x) 

Up 

Cf 



--2 + X 

X 



1 — X -\ — (5(1 — x) 



C, 



1+ x^ 



a;^ + (1 - xf 



(l-x) 



1 - X + -6{1 - x) 



(91) 



10 Appendix 2. Discretizations 

For the calculation of the weights used in the numerical analysis we follow closely ref. [6] 
also implemented in [7]. At each x- value, these authors use an approximation characterized 
by weights which are calculated analytically, together with an interpolation formula for the 
integration function. The method allows to monitor the singularity appearing in the small-x 



region and to achieve a very good numerical accuracy. The method speeds up in time the 
calculation by a large factor, but becomes tedious when moving to a higher order, since all the 
integrals have to be exactly discretized and the logarithms extracted in each sub-interval. Here 
we briefly illustrate the method as it applies to our case. 

We briefly recall the numerical strategy employed in this analysis. We define P{x) = xP{x) 
and A{x) = xA{x). We also define the convolution product 

J(x) . /' ^ (^-) P (^-) My). (92) 

J'^ y \y) \yj 

The integration interval in y at any fixed x- value is partitioned in an array of increasing points 
ordered from left to right (xo,a:i,X2, ...,a;„,a:„+i) with xq = x and Xn+i = 1 being the upper 
edge of the integration region. One constructs a rescaled array (x, a;/a;„, ...,x/x2^x/xi, 1). We 
define Si = x/xi, and Sn+i = x < Sn < Sn-i < ...si < sq = 1. We get 

^ f^«+i dy /x\ ^ /x 

i=o^-> y \y) \y 

At this point we introduce the linear interpolation 



7^J^^ y \y \y 



A{y) = 1 - JL^ A{x,) + y ^' A(x.+i) (94) 

and perform the integration on each subinterval with a change of variable y— > x/y and replace 
the integral J{x) with its discrete approximation Jn{x) to get 

Mx) = A{x,)-^ f '^P{y){y - s^) 
I- si Jsi y 



- J:A{x,)^^ r'^P{y)iy-s,_,). 

; 1 Si — 1 Si J Si V 



(95) 
Introducing the coefficients W{x,x) and W{xi,x), the integral is cast in the form 

n 

Jj^(x) = W{x, x)A{x) + J2 ^i^i^ x)A{xi) (96) 

where 
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(97) 



We recall that 
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(98) 



(99) 



as can can be shown quite straightforwardly. 
We also introduce the expressions 
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(100) 



Using the linear interpolation formula (94) we get the relation 



dyyAn{y) - xAn{x) 



X y 



y-x 



log(l - x)An{x) + An{x)InQ{x) 



N 



+ ^An{xi) {Jni{x) - Jnti{x)) 



(101) 



i=l 



which has been used for a fast and accurate numerical implementation of the recursion 
relations. 
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Figure 1: u, d, s and gluon distributions at Qi = 5 GeV 
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Figure 2: xu{x) and xd{x) evaluated with Qi = 2.0 GeV and Qf = 100 GeV with mx = 10 
GeV in the standard (reg) and susy evolution 
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Figure 3: Gluon distributions with Qi = 2.0 GeV and Qf 
GeV. The regular and the susy evolution are shown 
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Figure 4: xq^~^\x) (singlet) quark distribution evaluated with Qi = 5.0 GeV and Qf 
GeV with m\ = 10 GeV in the standard (non-susy) and susy evolution 
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Figure 5: Comparison of the supersymmetric gluon and gluino distributions for rrix = 30 GeV 
and Qf = 100 GeV. 
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Figure 6: Gluon distribution for 3 values of the gluino mass 100, 150 and 200 GeV and a final 
evolution scale Qf = 1 TeV. Shown is also the regular (reg) evolution of the same distribution 
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Figure 7: Gluino distribution for Q/ = 1, 2 and 5 TeV and mx = 100 GeV 
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8: Gluino and gluon distributions for very large final evolution scales Qf = 5 and 10 
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Figure 9: Gluino distributions for a varying rrix (40, 100, 200 and 250 GeV) with a fixed final 
evolution scale Qf = 1 TeV. 
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Figure 10: u, d quark distributions m\ = 250 GeV with a varying final evolution scale Qf = 5 
and 10 TeV. Shown is also the strange quark distribution with Qf = 5 TeV. 
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Figure 11: Dependence of the total 2-gluino cross section Cpp^ aa on the factorization scale in 
the QCD and SQCD cases {vri\ = 250 GeV). Shown are the factorization scales Qfact= 510, 
600 and 700 GeV. We show both the QCD and the SQCD resuhs 
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Figure 12: Dependence of the total 2-gluino cross section Ogg^w on the factorization scale in 
the QCD and SQCD cases (mx = 250 GeV). Shown are the factorization scales Qfact= 510 and 
700 GeV.We show both the QCD and the SQCD results. 
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Figure 13: Dependence of the total 2-gluino cross section Oqq^w on the factorization scale in 
the QCD and SQCD cases (m\ = 250 GeV). Shown are the factorization scales Qfact= 510 and 
700 GeV.We show both the QCD and the SQCD results. 



